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A SOLUTION OF ONE-DIMENSIONAL FREDHOLM INTEGRAL 


EQUATIONS OF THE SECOND KIND 
Ralph E. Gabrielsen 
Ames Research Center 


INTRODUCTION 


In the process of developing a rigorous numeric solution for the incom- 
pressible Navier-Stokes equations, the need has arisen for numerically solving 
Fredholm integral equations of the second kind in such a way that the complete 
relevance of the numbers obtained is known. Accordingly, this report addresses 
the problem for one dimension. The results are based upon the outstanding 
theory developed by Kantorovich. For one-dimensional Fredholm integral equa- 
tions it will be proven that the particular construction converges to the true 
solution. For purpose of clarity, it should be pointed out that Kantorovich 
(ref. 1) presents a related example for a periodic one-dimensional Fredholm 
integral equation; however, it is wrong in the sense that his proof does not 
comply with the necessary conditions of the theory applied. 


MAIN DEVELOPMENT FOR THE EQUATION 


For the equation 

f 1 

x(s) - X J h(s,t)x(t)dt - y(s) (1) 

0 

with h(s,t) continuous and periodic in both variables over the interval [0,1], 
and y also a continuous, periodic function over the interval [0,1], it will 
be demonstrated that an explicit system of equations of the form 

n 

x(t j ) - X £ AjhCtj.t^xUk) - y(tj) (2) 

converges to the solution of (1). In addition, it will be shown how the 
results can be directly extended to the nonperiodic case. 

Equation (1) will be considered as a functional equation in the space 
X = C of continuous, periodic functions on [0,1]. 

System (2) will be regarded as an approximate functional equation in the 
space X = R^, 

Kx = x - XHx = <py 


Also, equation (1) will typically be expressed in the form 

Kx = x - aHx = y 

Now consider the subspace ^ (or X) of X: Space of continuous functions 

on [0,1] that are linear on the partitions Aj_ where the T^'s (defined by 
» (i - l/2)/n, i » 1, . . . , n) are the end points of the interval Aj_. 

Define the mapping <J» 0 on X n onto X as follows: if x € X, then 
<j> 0 x - x, where x = (5^, • • Si = x(T^) , i - 1, . . ., n._ Since <J> 0 

defines a 1:1 mapping between X and 1^, exists. Norms on X and X are 

defined as follows: 

if 


xe X , ||x|| =» max |x(t)| 

o<t<i 


if 


x € X = Rjj , ||*|| - max 1 | 

i=i • . .n 


where 


x— ^2* • • • ’ 5 n ) 


Consequently, 

IK oil - IK‘H - 1 

Define the mapping $ as follows: 


if 


x e x , <J>x = 4 > 0 x 


if 

x € X , <j>x = [x(t x ) , x(t 2 ) » • • • * x(t^) ] , t^ — ^ 

•• Ill’ll * .. 11**11 ■ 1 

II *11-1 

Hence, it readily follows that I is a linear extension of <J> 0 . 

Let u> s (<5) denote the modulus of continuity of the function h(s,t) rela 
tive to s : 

co s (5 ) * sup|h(s + a,t) - h(s,t)[ (0 < s, t < 1, |a| < 6) 


2 


Similarly, u) t (5) is defined. Note, of course, that all parameters are 
restricted to the domain of consideration [0,1]. 

Now, with these preliminaries out of the way, we will quote a theorem of 
Kantorovich (ref. 1), which can now be directly applied to the problem in 
question. Consequently, most of what follows entails demonstrating that the 
conditions of the Kantorovich Theorem are satisfied. Before stating the 
theorem of Kantorovich, we will briefly discuss items necessary for the under- 
standing of the Kantorovich Theorem. 

We seek a solution of Kantorovich’s first equation (Kl) 

Kx = x - XHx = y 

in the B. space X, and the approximate solution x of Kx = x - XHx - Py in 
X C X where P is a projection operation from X onto X. 

Let <J> 0 define an isomorphism between X and the space X; now assume 
there exists a linear extension <j> of the operation <j> 0 to the space X. 

Let $ =* $ 0 P, P the projection operation from X onto X. Therefore 

P - (fro 1 * 


Hence, in (K2) 


<J> 0 (x - XHx » 4>7 lc t>y)3c » 4> 0 x 


yields the associated equations (K3) 

x - XHx = <(>y 


where 

H - ^H*; 1 , in X 

1. As shown by Kantorovich (ref. 1), if the following conditions can be 
satisfied 

A. for every x € X, ||h$ 0 x - <|>Hx|| < n||x|| 

B. for every x € X, there exists x S X such that ||Hx - xj| < ni||x|| 

C. an element y € X such that j;y - yj| < ^ 2 11^11 

and if the operator (I-H) has a linear inverse, and if 

q = |x|[n(l + |x|n 1 )||<j»; 1 || + n 1 ||<t»; 1 <DK||]l|K- 1 || < 1 
then (I-H)x = y has a solution x* for every right-hand side y€ X with 
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II. 


mils dh[ wi *»^ * - a ♦ WvlM IW1 ll K “ll 

Furthermore, if in addition, x* is the solution of (I-H)x = 

II x* - C*H S ?l|x*|| . 

and 

q ^ hjx| + n 2 ||i - h|| 

and 

P - (1 + e)|x|n||<j.; 1 K- 1 || + e[l + H+o^I- H) _1 4»(I - H) ||] 

III. Convergence results 
if 1. (I-H) has a linear inverse 

2. Conditions of Statement I are satisfied 

3. A, B, and C of Statement I are satisfied for n = 1, 2 

where 

lim nll^o 1 !!. =* lim n x || <J>7 1 $11 = l $ll “ 0 

n-*» n-*» 

then the approximate equation (K3) is soluble for sufficiently large 
the sequence of approximate solutions converge to the exact solution. 

II X* - x*|| < Qnllpo 1 !! + Q 1 n 1 ||$o 1 $|| + Q 2 n 2 || 'frazil 

where x£ = ^o^n* and Qx» Q» and Q 2 are constants. First, seek to 

||*Hx - h$ 0 x|| < n||s|| 


Case 1: Suppose 


Hx <= X 


Therefore, 


4>Hx ■ <J> 0 Hx 


Since 

H = 

<()Hx - H^gX = <J> 0 Hx - <() 0 Hx 


y , then 


n and 
Also , 


show 
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Therefore, 

||4>H5c - H<j> 0 x|| < || H - H|| ||x|| - n||*|| , for n = ||H 


Case 2: 


4>Hx ■ $ 


(Hx) <£ X , 
,t)x(s)ds , 


J h(s,t)5 
f h(s ,t 1 )x(s)ds , f h(s,t 2 )3 


)x(s)ds , . 


Seek to show that for every x € X , 

n 


||H4> 0 x - 4>H5cjl < n||x 


where 


H = d> 0 H4>7 X 


Hx 


j r 


h(s,t)x(t)dt; 


/ h < t n. 


4>Hx = J h(t 1 ,t)x(t)dt, J h(t 2 ,t)x(t)dt , 

'a o « 

<t> 0 x = [x(t x ) ,x(t 2 ) , . . x(t n )] 

H<J> 0 x = ^ [h(t 2 ,t 1 )x(t 1 ) + h(t 1 ,t 2 )x(t 2 ) + . . . + h(t 2 ,t n )x(t n ) , 

h (c 2 > 1 2 )x ( 1 2 ) + ti ( t 2 y 1 2 )x ( 1 2 ) ^ * 

h(tjj,t2.)x(t 2 ) + h(t n , c 2 )x(c 2 ) + . • . + h( , t^) x ( t^) ] 

.1 


t)x(t)dt 


-H<j> 0 x + 4>Hx 


h(t., ,t)x(t)dt £ h(t 1 ,t.)x(t ), J h(t 2 ,t)x(t)dt 
[ja 1 n i-x 1 1 *o 


i n r 1 

— 2^ h(t i » t . )x( t . ) , • • •» I 

i= i o 


h(t n ,t)x(t)dt 


* s E h(t n ,t^)x(t^) 


1=1 
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Let [0,1] be partitioned by the T^'s, i - 1, . . . , n where 
T ± = [i - (l/2)]/n, i - 1, .... n, with the length [T i . I .T^] =* | A ± | . Therefore 
there are n - 1 intervals of length A, with the end intervals = A/2. 


Let z(t) = h(tj,t). Consider 

r 1 i n 

J z(t)x(t)dt - - £ z(t i )x(t i ) 


n _T 


-Ef 


K+l 


|K-x % 


f ' A K+i 

[z(t K )x(t) - zCtjJxCtj^) ]dt 


K 


S u)A K (z)NI 

where 


E \ [ 2 (t K )x(t K ) -z(t K )x(t K )] 


K»i 


n ]-> z(t K )x( \ ) “ z(t K )x(t K ) 


T + T 
K K+i 

2 
n 


Z<t K* 

K-i 


X(t K+l ) + x(t K> 


- x(t R ) 


z<t K> 


K=i 


x(t K+i ) " x(t K } 


1 

Z ( t ! ) 

n 

2 

> 


z(t 2 ) 


2(t n ) 


x(t 2 ) ] + • • • + — 2 — [x(ti) - x(t Q ) ] 

[z(tn ) . 2(Ci)] [ 2 ( tl ) 


- z(t 2 )] + . 


x(t n ) 

+ 2 t z ^ c n-i^ “ z ^ t n^ 


< 1WL [ nUAK(2) ] < M UAk ( 2) 
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Therefore 


Hence , 


where 


1 

j z (t)x(t)dt 
o 



z(t K )x(t K ) 


< | co a (A K )|l»|| 


H H<j) 0 x - < n||.*|| 


n » j w s (a) 

For the nonperiodic case, an additional term comes in, due to the ends of 
the interval: 

C n 

J z(t)it(t)dt - 7 z(t K .)R(t Tr ) 

0 K*i 

Hence , 

l|H* 0 at - ^ hx H < n||*|| 

where n = (3/2)u s (A) + A||h|| for the nonperiodic case. 

Seek to show condition B. For every x € X there exists x € X such 

that 

II k* - *11 s "ilMI 


< [§ + ^llhlljw 


To this end, we will first show 


||z - 2 1| < a) 


Let 


V s * T j + i 

z (s) - 2 (s) = z(s) - n[(T J+1 - s)z(t j ) + (s - T.. ) z(t j+1 ) ] 


Therefore, 

z (s ) - 2 (s) = n(T - s)[z(s) - z(tj)] + n(s - T^)[z(s) - z(t j+1 )] 
| z (s) - 2 (s) | < n(T j+1 - s)oj 0 A^ + n(s - Tj)oj^ < A^ 
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Let 


2 * Hx 

z(s) - z(s’)| = J [h(s,t) - h(s’ ,t)]x(t)dt 


h(s' ,t) ]x(t)dt 


o 


* OJ s (y) ll X ll * I s' - s I < <s 


Therefore, 




Hence 


|| Hx - 2|| < 0 ), 


's 


(f ^||x||, where 


z = f^Hx 


For condition B, the proof for the nonperiodic case is identical. There- 
fore, condition B is satisfied with * a> s [(3/2)A]. 

To show condition C: Seek to show an element y S X exists, such that 


Therefore, the result directly follows by letting n 2 = (3A/2) ]/||y||. Simi- 

larly, this holds true for the nonperiodic case. Therefore, relative to 
page 7, conditions A, B, and C have been established for arbitrary n > 1, and 
it remains to show 

1. (I - XH) has a linear inverse 


Since n = (3/2)ui g (A) , - u g [ (3/2) A] , (n - (3/2)u> s (A) + |a| ||h|| , for the non- 

periodic case), and n 2 * <d[ (3/2) A] /||y|| , they all tend to zero as n -*■ ®; 
aiso, || (f.- 1 !! - ||*|| = 1. 

Hence, when X is a noncharacteristic value of the operation (I - XH) , 
then for sufficiently large n, the system of equations (k3) is soluble (with 
solution x*), and the sequence *o lx n converges to the exact solution of (1) 


II y - yil - n 2l|y|l 


From the preceeding section we see that 



2. q * | X | [ n (l + |x|n 1 )||*; 1 || + n 1 ||*; 1 *K||]||K- 1 || < 1 

3. lim nil*” 1 !! = lim n 1 ||4>“ 1 <)>|| = lim n 2 || <t>o ^ll = 0 
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CONCLUDING REMARKS 


Therefore, based on the Kantorovich theory, the approximate solutions 
Xq * s of the equations 

*n - Xfi n*n " ♦ o^ 

have been shown to converge to the solution x* of equation (1): 

x - XHx » y 


that is. 


||«* - x*|| -*■ 0 as n -*■ ® 
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